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Abstract 

We investigate a non-equilibrium reaction-diffusion model and equivalent ferromagnetic 
spin 1/2 XY spin chain with alternating coupling constant. The exact energy spectrum and 
the n-point hole correlations are considered with the help of the Jordan- Wigner fermion- 
ization and the inter-particle distribution function method. Although the Hamiltonian has 
no explicit translational symmetry, the translational invariance is recovered after long time 
due to the diffusion. We see the scaling relations for the concentration and the two-point 
function in finite size analysis. 
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The study of the systems with reaction and diffusion has been an attracting problem. Es- 
pecially, methods to determine the correlation functions and to observe the approach to the 
thermal equilibrium in one-dimensional models have been much developed recently. One of the 
recent success in this field is to obtain the exact correlation functions including the coagulation 
term by the so-called inter-particle distribution function (IPDF) method [|l|, ^ • If we appropri- 
ately tune the coupling constants, (i.e. the reaction-diffusion rates) in this approach, the master 
equation can be solved exactly. 

On the other hand, the (partially) antiferromagnetic XY spin chain with alternating coupling 
constant has been investigated as a toy model to consider the spin-Peierls phase transition or 
the Haldane gap. The thermal equilibrium model with the alternation has been studied by 
some authors ^. They have obtained the energy spectrum with the mass gap proportional 
to the strength of the alternation, the dimer correlation functions and so on. It has been also 
shown that the dimerizing process lowers the ground state energy. From those results, they 
have concluded that the uniform antiferromagnetic Xy-chain is unstable with respect to the 
distortion. 

In this paper, we investigate a thermal non-equilibrium model with alternating coupling and 
diffusion. It can be shown that the stochastic model and the ferromagnetic XY chain with the 
alternation are mapped on each other by linear transformations |2|, |5|. The basic tools to study 
these systems in this paper are the Jordan- Wigner method and the IPDF method. By the former 
method, which can be used in both the ferromagnetic and (partially) antiferromagnetic regimes, 
we can diagonalize the Hamiltonian exactly and calculate the energy spectrum. By means of 
the latter, which is available only in the ferromagnetic regime of the model, we obtain the exact 
forms of the energy spectrum and the correlation functions. Of course, the overlapped results 
obtained in the ferromagnetic region by both methods are the same. Since in the ferromagnetic 
case, the ground state occurs when all the spins are aligned, the time evolution of the correlation 
functions is our main interest in this paper. 

Hereafter, we investigate a spin-1/2 Xy-model in a magnetic field along the z-axis with 
alternation whose Hamiltonian is defined with the alternation parameter A and the diffusion 
parameter rj (> 1) by 



Hi,i+i = -V (??sf + ifif sf+i + sf + sf+i -V-^y (1) 

We assume that the length L of the spin-chain is an even integer and that periodic boundary 
conditions {sl+i = si) are imposed. 

Let us consider the exact diagonalization of the Hamiltonian (|I|) by the Jordan- Wigner 
transformation |§, ^, 0. It should be mentioned that this method is also available in both 
the antiferromagnetic (A < — 1) and ferro-antiferromagnetic coexisting (A > 1) cases as well as 
in the ferromagnetic (|A| < 1) one. Due to alternation, the unique band in the uniform model 
is split into two separated bands with an energy gap. Therefore, we introduce two kinds of 
fermions 62n-i and c„ corresponding to spins on the odd and even sites respectively. The spin 
operators can be represented by these fermions as 

(n-l n-l \ 

(n n—1 \ 

4j-1^2i-l +i-^Yl 4jC2j L 52n = (s2n)^ (2) 
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where the integer n runs from 1 to L/2. In terms of the fermions b and c, the boundary conditions 

are given by bi+i = ±&j and cl+i = ±Ci depending on whether M = I]^/i(4j_i^2j-i + 4j'^2j) 
is odd or even. Performing the Fourier transformation 



'J2n-1 



^ k ^ ^ k 



2 ■r ^ 27Tik r, 



and the Bogohubov transformation 

(bk,b'',Ck,c'') = (bk,P_k,Ck,c[^i,)'^k, 
where Sk is a certain (4 x 4)-matrix, we can rewrite the Hamiltonian (|^) as 

k 

Hk + Hl_, = a1+) {blbk - m'^) + aI") (4c, - c'^tc'^) . 
The energy eigenvalues A^^^ are given by 



\ 



cos2 f ^^"j + A2 



sm 



+ 



r] — 1] 



where the value of momentum k depends on M as 



k 
k 



L 



1, for7\A = odd, 

2 



0,1,2,--- 

13 LI 
2' 2'"'' 2" ~ 2' 

Therefore, the ground state is given by the half-filled state in this basis like 



for M 



even. 



mi] m 



\vac) 



with bk\vac) = b''\vac) 



Ck\vac) 



c \vac) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



and the energy spectrum has a gap 2r]A{r] + r]~^) at 2'Kk/L ~ 7r/2 provided A 7^ as is shown 
for that without the diffusion (i.e. the magnetic field) Contrary to the antiferromagnetic 

case, the energy gap does not appear near the Fermi surface in the ferromagnetic case. It follows 
that the ground state energy remains a monotonically increasing function of the alternation 
parameter A even if we add the elastic energy of distortion to the Hamiltonian. It means that 
the instability due to the dimerization cannot be observed in the ferromagnetic case. 

On the other hand, as is shown in [|2|, ^, the spin Hamiltonian can be mapped onto a 
stochastic model including coagulation and decoagulation processes by regarding the spin-up and 
spin-down state as sites occupied by a particle (denoted A) and empty (denoted (/>), respectively. 
If the i-th site is occupied (empty), we label Ui = 1 (0). The stochastic model to be considered 
hereafter is defined by the transition rates Wa/sifJ^,!^) (a, P, fJ-, i' £ Z2) and arbitrary positive 
constants a±. Only the following processes are allowed. 

1. Diffusion 

A + <j)<^4) + A at the rate a±wii{01) = a±wii{W). 
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2. Coagulation 

A + A —>■ A + 4) and A + A^(j)-\-Aa,t the same rate a±woi{W) = a±wio{Ol). 

3. Decoagulation 

A + (p ^ A + A and cp + A^A + A&t the same rate o-i-tt;oi(ll) = a-i-iz;io(ll)- 

The above rates should be interpreted as below. For example, let us assume that the i-th and 
[i + l)-th sites are occupied and empty, respectively, and that i is an even (odd) integer. After 
an infinitesimal time interval dt, the sites will be occupied by (p and A, respectively, at the rate 
a^wiQ{01)dt {a-Wio{01)dt). We restrict the rates further as 

w;ii(10) = wn{Ol) = woi{W) = u;io(01) = 1, 

z/;oi(ll) = u;io(ll) =r/2-l, (9) 

which makes it possible to solve the model exactly by the IPDF method as will be shown later. 
Because we can fix the normalization of a+ and a_ as a+ + a_ = 1 by rescaling time, we put 
a± = (1 lb A)/2 with the same A as in (|^). The basic quantity in the stochastic model should 
be the probability distribution function P{g_;t), which is the probability to find the system in 
the configuration a_ = (fii, ■ ■ ■ ,(Tl) at time t. Defining the Hamiltonian by 

L 

H = ^'^a(^_y Hjj-^-i with a(_-)odd = a_, a(„)evon = a_|_, (10) 



with periodic boundary conditions i = i + L, we can write the master equation of the probability 
distribution function in terms of a Schrodinger equation with imaginary time 



i=lPj,Pj+i=l,0 ' 
xP(fJi,- • • ,(Tj_i,Pj,/)j+i,crj+2,- • ■ ,crL]t). (12) 



The Hamiltonian above given can be transformed into that of the spin system (|l]) by the same 
mapping proceeded by (2|, |5[. Performing a similarity transformation Hi = U~^HU with the 
matrix 

..(v^;).....(^ J), as) 

and a rotational transformation H = R^^HiR with 

R = e^p{ies\) • • • exp(i6'4) (tan = rf - 1) , (14) 

we obtain the spin Hamiltonian (2). 

To calculate the correlation function from the master equation ([T^), we use the IPDF 
method g |, § -0]. For simphcity, we assume a Gaussian initial condition for general N 
and independent of M as 

^(ffl * = 0)^'^M,0 • • • 5<TM+N-lfl = P^, (15) 
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where p (0 < p < 1) is the probabihty of a site to be empty. If we define the "hole length 
probabihty" (HLP) 0,{m;t), which is the probabihty to find a string of empty sites with the 
length m at time t, by 



n{2n + l;t) = ;^P(o:;t)(5<,„,o---<5, 



J^W(2n;t) = 5]P(ffit)<5,„,o---<5.„^2„_i,o, (M = even/odd), (16) 



we can easily verify that the above fi's do not depend on the starting site M explicitly because 
of the translational invariant initial condition (|l5|). Note that there are two kinds of the HLP's 
for even hole-lengths, depending on whether the starting point M is even or odd. With 
these preparations, we can rewrite the master equation simply as 

d 

n{2n-l;t) = if- 



dt 



_J](+)(2n; t) + a^n^~\2n; t) - (1 + ??^)0(2n - 1) 



a+9.^^\2n-2;t)+a^9.^-\2n-2-t) for 1 < n < ^, 



d r 

— J7W(2n;t) = 2a:F [ry20(2n-l;t)-(l + r?^)OW(2n;t) + 17(2n + l;t) 



L 



for 1 < n < 1. 

- - 2 

The solution of Eq. ( [T7| ) is given by 

L /2— 1 

f](2n-l;t) = ^4^^??-(2"-i)e^^'^'*sinf^^^^^^A + 

1=1 s=± \ L J 

+i32„-ie-('''+^)* + V(2n-l), 

QW(2n;t) = F + + 1) ^,(^),sin (2vrn//L) 

^ ' ^ ^ I I cos(7r//L) 



(17) 



+ 



1 — a± 



1 - 2a± r?2 + 1 
i(±) 



V(2n). 



(18) 



In the above expression, Ai' and 52„ _i are complicated coefficients depending on L, A, ry 
and / and the zero-mode function is 



The energy spectrum X^^^ takes the same form as that in (^). 

With the above solution, we can investigate the finite size scaling of physical quantities. 
From now on, we set r] = 1, i.e. the massless regime. First, we consider the finite size correction 
of the concentration c{t) = 1 — 0,(1] t) H]. It is not difficult to ensure that we can fix a parameter 
z = 2(1 — A^)t/L^ finite in the scaling limit L ^ oo and t — > oo. After performing a modular 
transformation z — > —1/z, the asymptotic form of the concentration reads 



1 



1 — r] 



-2L 



(1 - p^)r]~^"' +p^- r/-2^ 



(19) 



c{t) 



1 



27r(l - A2)t 



32(1 - A2)f 



l + 6p + p^ 



+ A 



3-2p + 3p^ 
(l+p)2 



+ 0(t 



-5/2] 



(20) 
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The finite size analysis for the two hole probability Q^^\2;t), which is identified with 
the probability to find a dimer on a link for a±, can be done similarly. Due to alterna- 
tion, r2(+)(2;t) and ^l^~\2;t) do not coincide for general t in spite of the initial condition, 
O(+)(2;0) = n^-\2;0) = However, because of the diffusion and coagulation, we will observe 
il(+)(2;oo) = n^-\2;oo) = 1, which indicates the recovery of the translational invariance. For 
example, let us consider the ratio of the two probabilities i}^^\2;t) and i}^~\2;t). After a 
lengthy calculation, its finite size scaling is shown to be 

As we have seen, both the Jordan- Wigner method and the IPDF approach are available in 
the XY model with diffusion and alternation. An energy gap proportional to the strength of 
the alternation is observed. We have also seen that the translational invariance is recovered 
by the diffusion beginning with a translational invariant initial condition. This recovery, e.g. 
0(-)(2;t)/fi(+)(2;t), is scaled by the same variable z as that of the concentration. We have 
obtained the exponents of the concentration and the two hole probabilities in the scaling limit 
L, t — > oo. 

We can imagine some extension in this field. First, the coupling constants do not have to 
change on every other link in our approaches. That is, it is possible to change the coupling 
constant on every third (fourth and so on) site without loss of solvability. In particular, it would 
be interesting to consider the relation to the model with randomly changing coupling constants 
|lT|| . Secondly, contrary to the IPDF method, the Jordan- Wigner method can be also applied for 
the (partially) antiferromagnetic alternating model. Although we did not succeed in finding the 
corresponding initial condition in the Jordan- Wigner approach so far, we expect that a parallel 
discussion can be done. 
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